Hidden order in bosonic gases confined in one dimensional optical lattices 
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We analyze the effective Hamiltonian arising from a suitable power series expansion of the over- 
lap integrals of Wannier functions for confined bosonic atoms in a Id optical lattice. For certain 
constraints between the coupling constants, we construct an explicit relation between such an ef- 
fective bosonic Hamiltonian and the integrable spin-S anisotropic Heisenberg model. Therefore the 
former results to be integrable by construction. The field theory is governed by an anisotropic non 
linear a-model with singlet and triplet massive excitations; such a result holds also in the generic 
non-integrable cases. The criticality of the bosonic system is investigated. The schematic phase 
diagram is drawn. Our study is shedding light on the hidden symmetry of the Haldane type for one 
dimensional bosons. 
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I. INTRODUCTION 

Systems of cold atoms trapped in optical lattices pro- 
vide a remarkable tool to simulate quantum many body 
physics in engineered quantum systems 0, HJ. In this 
context new perspectives are provided by the possibility 
to handle atoms or molecules with large dipole moment 
0. For these systems, the Bose- Hubbard models can 
be employed provided that processes arising from the 
long range interaction are considered. The correspond- 
ing 'extended' Bose-Hubbard Hamiltonian can display a 
variety of quantum phase transitions between supcrfluid 
states to Mott Insulator (MI) ones (f.i. see Q). 

Recent investigations by Altman and coworkers indi- 
cate that a strongly interacting bosonic system can dis- 
play a further gapped phase with 'exotic' order [f|. Based 
on numerical analysis of a Bose-Hubbard type model, 
such phase was demonstrated to be characterized by a 
non local hidden order of the Haldane type. Originally 
discovered in spin systems @, [3 , Haldane phases are be- 
lieved to play a crucial role in many different contexts 
including QCD and High— Tc superconductivity pd| . 
Recently the hidden order was noticed also in different 
types of electronic insulators fl2l. [l3j . 

The Haldane order in the bosonic systems was evi- 
denced by rephrasing the physical concepts provided to 
study of the hidden order for spin systems. Although 
Altman and coworkers constructed a variational scenario 
in order to interpret their results, it would be desirable 
to trace how the bosonic hidden order emerges from the 
Haldane order in spin systems. This is a major purpose 
of the present work. 

Based on exact methods we construct a precise relation 
between the bosonic and spin hidden orders. Relying on 
it the quantum criticality of the bosonic system is inves- 
tigated. We construct the non-lincar-cr model (NLcrM) 
describing the system at long wave lengths. The phase 



diagram is investigated within the saddle point approxi- 
mation of the NLcrM. Our results are shedding light on 
the structure of the hidden order in bosonic systems, sug- 
gesting that the bosonic excitations in the Haldane Insu- 
lator (HI) arc of triplet nature. 

The paper is organized as follows. In Sec. II we derive 
the model Hamiltonian from the microscopic dynamics 
of a dilute system of interacting dipolar bosons. The 
emerging Hamiltonian defines a type of extended Bose- 
Hubbard model of interacting bosons in the lowest Bloch 
band. In section III, we establish the mapping between 
the bosonic model and the integrable higher spin-S* XXZ 
Hamiltonian [li] (see Eqs.(l4l)- (flTJ|) ). Because of the inte- 
grability of the spin model, the bosonic Hamiltonian is in- 
tegrable by construction as long the microscopic parame- 
ters entering into bosonic Hamiltonian fulfill certain rela- 
tions (see Eq. fTTj) ). Relying on this result, the quantum 
criticality of the bosonic system is investigated in the Sec. 
IV where we construct the non-linear-cr model (NLcrM) 
describing the system at long wave lengths. The phase 
diagram is constructed, beyond integrability, within the 
saddle point approximation of the NLcrM (see Fig.([T]); we 
note that the results are not restricted by Eq. lfTTj) ). In 
section IV we discuss and suggest experimental protocols 
to detect the Haldane insulating phase. Finally we draw 
our conclusions. 



II. THE BOSONIC MODEL HAMILTONIAN 

We start by analyzing the order of magnitude of the 
amplitudes involved in the Id lattice Hamiltonian for 
trapped bosonic atoms; we include dipolc-dipolc inter- 
action. We will obtain an effective model with density- 
density interaction and higher order hopping processes. 

The general Hamiltonian for the atomic (of mass m) 
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gas reads H = Hq + Hi nt . with 



^V 2 - V(r) 
2m y ' 



*(r), (1) 



H mt = J drdr^ + {r)i! + {r)V lnt .{r-r)i>{r )i!{r) , 

where V results from a combination of harmonic confine- 
ment with the optical lattice: V = Vharm + Viatt- We con- 
sider a 'cigar-shaped' configuration Vharm — muj 2 (x 2 + 
7 2 2/ 2 + 7 2 z 2 )/2 where the harmonic potential has a fre- 
quency u> along the x direction and is much more confined 
of an anisotropy factor 7 3> 1 in the y, z directions. The 
Id lattice, Viatt = sE r sin 2 (nx/a), is arranged along x. 
E r = (for) 2 /2a 2 m is the photon recoil energy, a is the 
lattice spacing and s measures the optical lattice depth 
in terms of E r . 

The interaction potential Vint. = V sr + Vdd contains 
two terms: an 'onsite' short range potential V sr (r) = 
4irh 2 aBBS(r)/m characterized by the s-wave scattering 
length olbb\ and a 'long range' anisotropic dipole-dipole 
potential V dd {r-r) = ^^ 2 {l - 3 cos 2 0)/(47r|f- r | 3 ), 
/i being the atomic magnetic dipole (ju is the vacuum 
magnetic permeability) , and 9 being the angle of r — r 
with the dipoles orientation. 

The bosonic field operator can be realized through 
Wannier functions w(f — f\) localized around ff. 
i5f(r) = Y). bjwjr — rj), where fej annihilates a boson at 
the lattice site i. In this formalism the Hamiltonian 
reads H = V, , U^bj + J2 i]kl tij. k ib\b\b k b h where Uj 
and Uj-ki are the standard Wannier functions integrals 
They can be expanded as power series of the so called 
'lattice attenuation parameter' e = exp(— a 2 /Al 2 pt ) = 
cxp[-(7r/2) 2 x /i], with l opt = a/ins 1 / 4 ) < a 0; w{r-n) 
are assumed to factorize: w(r — f$) = w(x — Xi)w(y)w(z) 

and approximated as w(v) = (■n~ 1 / A l„ 1 ^ 2 )exp(— v 2 /2l 2 ), 
v = x 7 y,z, with l x = l opt and l y = l z = Zj_ = 
(h/m-fLu) 1 / 2 . 

The effective lattice (grand canonical) Hamiltonian 
emerging from the analysis above is H = E r H{, with 

H b = -K^rn -ty]bjbi + i + Uo^niirii - 1) 

i i i 

+U 1 '^2nin t+1 - t e y](n, + n t+1 )b\b t+1 + 

i i 

^$>I) 2 (&i+i) 2 + --. , (2) 



where rii = b\bi is the number operator. Terms in the 
ellipses involve higher powers in e (see also [3] ) . The first 
two terms, with k being the chemical potential and t the 
nearest neighbor hopping (in units of E r ), are the only 
low order terms from H since the matrix elements 
decrease as el 4- - 7 '! 2 . U = (2/7r) 3 / 2 a BB as 1 / 4 // 2 ^ is the on- 
site interaction neglecting the renormalization due to Vdd 
(that is of order e 2 ). Defining Idd = mfion 2 / (2air 3 H 2 ) the 



contributions of Vdd to the integral i«+i;n+ij the coupling 
constants in ([2|) read 



to = C/ e 3/2 



2t p = {U + Idd)e 2 



(3) 



We comment that the (form of the) second quantized 
Hamiltonian |(5J) is not affected by the gaussian approxi- 
mation of the Wannier functions which might only mod- 
ify the expression of the coupling constants in l[3"]). 

Besides the density-density interaction U± , provides 
for correlated t c and bosonic pair t p nearest-neighbors 
hopping processes. It is important to note that exper- 
imentally there are different ways of manipulating the 
relative strength of the coupling constants e. g. it is pos- 
sible to change the optical lattice parameters [15[ , or use 
appropriate Feshbach resonances [lj| to separately tune 
Idd and Uq, independently from the expansion in terms 
of the parameter e. This allows for the possibility of ex- 
ploring a large portion of phase space. 

Below we demonstrate that the Id bosonic Hamilto- 
nian (Eq. ) is obtained as large S limit of the intcgrable 
spin model generalizing the XXZ model to higher spin 
PjJ ■ This will prove that the model ^ is integrable (see 
however the discussion below HT 



III. 



THE INTEGRABLE HEISENBERG MODEL 
FOR HIGHER SPIN 



Integrable XXZ models for higher spin [l9| were 
intensively studied in the framework of the quantum in- 
verse scattering method [2(|. The Hamiltonian (E r is 
assumed as the energy's untit) is obtained as logarithmic 
derivative of the transfer matrix 21 1 



H 



XXZ 



(4) 



where 



ip a [x] = T' a (x)/T a (x) 

00 

r a (x) = (l-O^IIl 1 



n=0 



(5) 



reducing to the ordinary gamma function for a — > 0. The 
quantity J{a) is related to the (coproduct A of) Casimir 
operator C5 = S + S~ + svahaS z smha{S z + l)/sinh 2 a 
of the quantum algebra su a (2) underlying the integrabil- 
ity of the theory: 



sinhaJ(a) sinh a(J(a) + 1) 
2 sinh 2 a 



(6) 



We remark that the anisotropy a enters into |?|) deform- 
ing both the gamma function and the 'representation' 
J (a) of su a (2). 
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The limit of large 5, with small aS is interesting 
for us. The first 'effect' of such a limit is that the 
ipa(x) m © reduces to the ordinary digamma function 
ip{x) = r (a;)/r(a;). In such a limit the quantity J(a) 
can be obtained explicitely by resorting to the explicit 
expression of ACs 



AC 



(M+i) 



n aS z 



i+1 + 



. cosh aS cosh a(S + 1) 

smhai,- smnai,",, = 

smn a 

smhaS sinha(S' + 1) 



■ cosh aS* cosh aSf +i 



sinh 2 a 



-aS z 



(7) 



By retaining the terms up to the second order in aS we 
obtain 



AC* 



(*>*+!) 
S 



_ QX QX | QV QV i \QZQZ 

+ y((Sf) 2 + (ST +1 ) 2 ) + S(S + 1), (8) 



where A = l + a 2 fcwithfc = 5(5 + 1) + 1/2. By exploiting 
its relation with ACs, J{ a ) can be written as 1/5 per- 
turbative expansion: J(a)i.;+i = A{a)i_i+i + B(a)i^ + i 5 
to finally give 



J{a) 



i,i+l 
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+l/2[-3a 2 /8 + (1 - 3a 2 /8) (b\b l+ i + h.c 
-(l + 9a 2 /8)(h l + n l+1 )} (9) 

where we used the Holstcin-Primakoff realization of su(2) 
[22l |. The bosonic model (J2]) is obtained first employing 
the Stirling formula for the asymptotics of the digamma 
function ip(x) ~ hi(x), then by expanding the equation 
resulting from (J4|) at second order in 1/5: 



Hxxz = H b + const. 



(10) 



where const. = iVln[5(2 + 3a 2 /4)] - Nip(l). Transla- 
tional invariance has been assumed. The spin 5, and the 
anisotropy a are obtained by comparison of (|10p with 
©: 5 = t/{2U + Ui), a 2 = 2(2U - f7i)/(3E/i), with 
the constraints 



. _Ui 2t c - U 

c ~ 2 ' p ~ 2 
(2[/ + [/ 1 ) 2 



8(2C/ - 3E7i) 



8/- 



(St-Wf/o + H^i) 2 



U 2 Ux = (2U + U!Y 



(11) 



We note that only one parameter results to be adjustable 
in @ achieved from ; we remark that the restrictions 
(|TT|) can be achieved by tuning the relative strenght of 
contact versus dipole interactions. The resulting one- 
parameter- Hf, is integrable by construction (see also [23|); 
the exact solution will be studied elsewhere. In the 



isotropic case a = 0, is obtained from the Faddeev- 
Takhadjan-Tarasov-Babujian model (24[. In this case: 
Ui = Uo = t 2 = 2t c = At p H3. 

The ground state of ^ is a singlet: 5f ot = 0. For imag- 
inary a the spectrum is gapless. For real a it was proved 
that the excitations are gapped [l!|. Given the relation 
pop , it is intriguing to study the low energy spectrum of 
Hxxz in the limit of large 5. This is what we are going 
to do by exploiting the NLcrM. 



IV. NON LINEAR SIGMA MODEL 

As customcrly, the large 5 limit is combined with the 
continuous limit a — > 0. The action S = iSu[Q] + 
J dTdxH.[Q](r, x) is obtained within the spin coherent 
states path integral formalism: Ti[Vl] = (Q,\H\{T) and 
iv[n] is the Berry phase (22I]. 

The semiclassical continuous Hamiltonian arising from 
gradient expansion of (J!]) is H = {Q(x, t)\ACs\Q{x, t)) + 
. . . where the ellipses indicate terms with higher powers 
of a and 1/5 (or combinations of thereof). TC has A — D 
form with D = a 2 ; the gradient expansion procedure 
spoils the integrability of the lattice model (ITIB. since the 
t c and t p terms turn out of higher order in q[26|: never- 
theless integrability manifests in the field theory in form 
of the restriction A = 1 + kD, with k = 5(5 — 1/2), aris- 
ing from the lattice theory. Generic values of A, D can be 
considered through the 1/5 expansion of the A — D lat- 
tice model (that is not integrable), providing a different 
paramctrization of the spin parameters in bosonic terms: 

4i+ 1 

5 = — — , A = E/i , D = Uo . (12) 
o 

The large scale behavior of the system is captured by the 
large 5 realization of the coherent state spin variables 
in terms of staggered magnetization nj(r) and quasi- 
homogeneous lj (t) fluctuating fields through the Haldane 
mapping 0] 



(n\S 3 (r)\n) = (-iyn J (r)\l 



\lj(r)\ 2 , h{r) 



S 2 



s 



, (13) 



where ^(r) 2 = 1 and fj(r) • n,-(r) = 0. The NLcrM de- 
scription of (J3|) is obtained following a variation of the 
procedure originally adopted in 0, [29|, [3l[ to deal with 
the 5 = 1, A — D model. Due to the anisotropy of the 
A — D model, we separate the staggered magnetization 
and its fluctuation in their perpendicular and transversal 
components: rij = (n±,n z ) and lj = {l±,l z ). In order 
to obtain the Lagrangian of the NLcrM one needs to in- 
tegrate out the I field in the Hamiltonian together with 
the terms coming from the Berry phase. After taking the 
continuum limit one finds 



5 2 |cW| 2 + ^|d T nJ 2 



(14) 



a" 
~2 



5 2 A (d x n z y 



MS z n 



2„2 
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FIG. 1: Saddle point phase diagram for the NLcrM. Based 
on the behavior of the gaps (inset), the different phases are 
identified through the indications we reported in the text; the 
small offsets in the decay of A z and Ax is artifact of the 
saddle point approximation. The plot corresponds to t — 
15/4. Along the orange continuos line U\ = 1 + kUo the 
lattice model ([2} is integrable. The green dotted-continuos 
line indicates the degeneration of the gaps: A z = Ax (a 
similar phenomenon was evidenced for S = 10, [2SL l3ll]L 
Inset: Behavior of the gaps A z (dashed lines) and Ax for 
Uq = 1 and for different t (leftmost blue lines correspond to 
t = 15/4 while red lines to t — 23/4). For increasing t the HI 
is shrunk. 



with M = 1 + S(S - 1/2)U - Ux, plus the topological 
phase associated with the Neel field n = (n±,n z ) 



If 

T = — j drdx n ■ (d T n x d x n) . 



(15) 



T can assume only integer values for it is the winding 
number of the mapping n : M CO mp S> 2 which is classi- 
fied by the second homothopy group 7r 2 (S> 2 ) = 7L. T con- 
tributes in the partition function as e 2lT7rS , with T £ 7L, 
so that it does not contribute for integer spins. In Eq. fM]) 
the Neel field satisfies the constraint |n| 2 = 1. The coef- 
ficients are 



(ji- Mn 2 z )(2 - Mn 2 z ) 
u — 2n 2 (l - Ux) 
(fi - Mn 2 z ){2 - MnlY 



(16) 



with fi = 1 + Ux + kUo. The coefficients provide an ad- 
ditional interaction between the components of the Neel 
field. For M = 0, Eq. Q3|) is the sum of an 0(2) NLcrM 
in the field nx and a scalar model in n z . The former 
describes a free Gaussian model with a bosonic compact- 
ified field which is known to be a conformal field theory 
with central charge c = 1; the latter is also integrable. 
One of the most striking features of the Haldane phase is 
that the gap is provided by a triplet structure (7l.l27l|. The 



masses of the particles Aj_ and A z belong to the sectors 
S z = and S z = ±1 respectively. For the case 5=1 
A z and Ax play the role of the excitations of a confor- 
mal field theory with c = 1/2 and c = 1 respectively [29|. 
Eventhough the Haldane order in the bosonic systems 
was evidenced through a non local 'order parameter' de- 
fined in analogy with the spin string order parameter]^], 
the phase diagram can be studied also by analysing the 
low energy properties of the system and specifically the 
particle-hole S£ c and the neutral S£ n energy gaps@. We 
shall see that A z and Ax play the role of S£ c and S£ n 
respectively. 

To obtain such quantities from the continous field the- 
ory we treat (|X4|) by saddle point approximation (S = 1). 
The method assumes that < n 2 z >:= £. This means that 
|| (Srij) 2 — SnjSrij + i\\ <C 1 with 6n = n — fi, indicating 
that only sparse charge-fluctuations around 'diluted DW 
states' (so called 'zero defects states' [3(|) are considered. 
The constraint on the field n is taken into account by 
introducing a uniform Lagrangian multiplier r\ 



C = — 



a 
2 

a 2 
2 



-|&nx| 2 

9± 

—A (a x n z ) 



|9 T nx| 



v±g± 

(d T n z y 
v z g z 



'2„2 



MS'n 



+ r)(x,T)(\n\ 



(17) 



where v±/g± = S 2 , v z /g z = S 2 \ and l/(w i c/ i ) = a/2 
with i = {_L,z}. There are two poles in the propaga- 
tors of the effective action, providing the expressions in 
terms of £ for the gaps in the z and in the perpendicu- 
lar direction, respectively A z and Ax, to be determined 
selfconsistentlv 3l|. In order to derive the expression for 
the gaps we need the calculate the propagator. We sub- 
stitute in the action the Fourier transform of the field n 
and of the Lagrange multiplier rj, respectively n and t). 
After performing a Gaussian integration on the field n, 
the effective action becomes: 

S[fj] = J2 {HK- 1 )"} (q,n;q',n')-ri{0,0) , (18) 

q,n,j 

where (K^ 1 )^ are the diagonal entries of the inverse 
propagator 

(K~ 1 )" = 



1 

JL 



Snn'S, 



2g jVj 



(ft 2 +q\)- -±{Aq,An)+n S, 3 



(19) 



with Aq := q - q' , An := n - n' and j = {1, 2, 3}. The 
notation Vj, gj stands for v±_ g± if j = {1,2} or v z and 
g z if j = 3. fl n are the Matsubara-Bose frequencies, 
f2„ = 2iTn/(3. (AT -1 ) 33 provides an expression for (. 

With the ansatz fj(q,n) = (/3L)S n oS q o rj, we can derive 
the saddle point equation (dS/drj = 0). For low tempera- 
tures ((3 — * oo ) and in the thermodynamic limit (L — > oo) 



5 



one finds 



A6 



+ c 



+ 



A 2 



i 



(20) 



where the ultraviolet cutoff A is a free parameter to be 
determined by fitting the gap to the known value for the 
isotropic point A isot r = 0.41048 where £_l = £ z . The 
notation stands for £~ 2 := — with 7 = {_L, z}, £~ 2 := 

^J 2 /A and ^~ 2 := 2A u /g ~ . The expressions for the gaps 



A z =V z y/^ T 



(21) 



We can solve numerically the self-consistency equations 
(f20f for different values of U\, Uo and of the spin S (or 
equivalently i). From the numerical study of the gaps 
one can draw a (qualitative) phase diagram (see Fig[T]). 

We see that A z and Aj_ display a behavior that is very 
similar to S£ c and S£ n respectively (see FigQ]). In such 
a view, the phenomenology of the bosonic excitations in 
the HI would arise from the triplet nature of A 2 and A^ . 
In particular, the line S£ c = S£ n that was evidenced nu- 
merically in Ref. Q can be interpreted as the degeneracy 
of the triplet excitations of the field theory (|T4|) (that 
includes the 0(3) model A = [13). 

Further insights on the criticality of the bosonic system 
can be obtained by adapting the results for the spin S, 
\ — D modelJH]. Specifically, the onset to the Hi-density 
wave phase is suggested to be of the Ising type (second 
order), c = 1/2 and it is characterized by A z = with 
Aj_ 7^ 0; the MI-HI phase transition is with c = 1 and it 
is caused by A z ^ with Aj_ = 0. Similar behavior of 
S£ c and S£ n was noticed in Q (except that also 5£ c = 
was displayed to vanish at the MI-HI transition). The 
density wave (DW) phase-MI transition is predicted to 
be of the first order; the c = 1/2 and c = 1 lines meet in a 
point that is governed by the continuos limit of the a = 
integrable theory, with c = 35/(1 + S), S = l/(2i)[Hj]. 
The phase diagram of the system is summarized in Fig [I] 



V. EXPERIMENTAL FEASIBILITY 

Although a detailed analysis would go far beyond the 
scope of the present paper, we would like to sketch a 
possible method that could serve for detecting the HI 
experimentally (complementing the observations in 0]). 
The basic idea relies on the detection of the atomic cur- 
rent tracing back to the early experimental evidence of 
MI phase J3J]: the MI conductivity is probed by apply- 
ing a washboard potential to the lattice. A resonance in 
the atomic conductivity appears when the tilt between 
adjacent sites reaches the energy gap Uq i-e. when it 
is resonant with the particle-hole pairs excitation energy. 
Because of the peculiar solitonic non local order [3(| of the 



HI the resonance peak in the conductivity is narrowed, 
as already noticed with a different experimental situation 
(parametric resonance) by @ . Placing the bosonic chain 
in a washboard potential, the distinct feature of the HI 
phase would be a dependence of the atomic current on 
the length of the chain (this is the analog of the diffusive 
spin transport evidenced in Haldane compounds [35j ) . We 
observe that, being an incompressible phase, the HI is ex- 
pected to be robust to the parabolic confinement when- 
ever the induced total energy offset (between center and 
the trap edges) is of the order of the energy gap. Fol- 
lowing the numerical indications provided in [f| such a 
gap is estimated to be of the same order of magnitude 
of the MI one: A ~ t. Therefore we can estimate the 
magnitude of the allowed harmonic confinement to be 
mu) 2 x 2 M < t = 2seE r , xm being the size of the conden- 
sate along the optical lattice. It would be interesting 
to trace the analogous for the HI of the "wedding cake" 
structure that appears in MI in presence of strong con- 
finement. 

Finally, ring shaped (with circumference L) optical lat- 
tices with twisted boundary conditions [36j could be ex- 
ploited to evidence the HI in an expansion experiment 
(we note that the energy offset due to harmonic confine- 
ment can be minimized in this case). In fact, adapt- 
ing the results obtained for the spin diffusion in twisted 
Heisenberg integer spin rings [37], the atomic density cur- 
rent would display a characteristic parametric depen- 
dence on the boundary twist: a sawtooth like behaviour 
for correlation lenght £ 3> L, and it would be exponen- 
tially suppressed, with sinusoidal oscillations for £ <C L. 
This could be a fingerprint of the Haldane gap for (finite) 
ultracold atomic systems. 



VI. CONCLUSIONS 

By exploiting a suitable expansion of the matrix ele- 
ments in terms of the lattice attenuation parameter e«1 
we derived an effective model for bosonic atoms in a Id 
lattice (J5J). Additional terms enter the Hamiltonian re- 
spect to the standard Bosc-Hubbard model. For certain 
choice of the coupling constants, the model results to 
be integrable through a mapping with the spm-S-XXZ 
integrable model (see Eq.([10P). We note that the possi- 
bility of independently tuning the onsite- interaction rela- 
tive to the density-density interaction, is available in the 
bosonic model Hamiltonian only in presence of dipole- 
dipolc interactions. This is precisely what allows degen- 
erate quantum gases to explore different portions of the 
phase diagram beyond the MI region. 

The direct relation between the spin and bosonic pic- 
tures are exploited to investigate the critical properties 
of the bosonic systems. From the present work, it is sug- 
gested that the correct context is provided by the spin 
S, X-D model. The HI is investigated by studying the 
continuos field theory arising from gradient expansion 
of the lattice model: the Lagrangian has the form of 
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a NLcrM with further interactions between the compo- 
nents of the Neel field (sec (fT4|) . (fTS]) ). We comment that 
the terms t p and t c entering the ^ do not contribute 
to the NLctM formulation indicating that such terms are 
irrelevant for the criticality we studied. Such an effect 
is a physical manifestation of the Affleck analysis [2(| . 
The difference between intcgrablc and non integrablc lat- 
tice theories is reflected in a different paramctrization of 
the coupling constants entering the NLcrM. Interestingly 
enough the integrable parameterization of the continous 
model (fT4|) we found should be tightly related to the 
so called pricipal chiral fields [38l| indicating an emergent 
SU(2) x SU(2) symmetry of the bosonic Haldane phase. 
We notice however that the phase diagram of the system 
(|14|) is investigated, beyond the integrability, by means 
of the saddle point approximation; the integrablc case 



might need separate discussion. The Haldane phase is 
found for any finite range of the interaction. This should 
be relevant for experiments where the scattering length 
can be tuned around zero thus evidencing the interaction 
between light-induced dipoles 
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